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We report on a search for the CP asymmetry in the singly Cabibbo-suppressed decays —> 
K~ K^-K^ and in the resonant decays ^vr^ and — > K*'^K'^ based on a data sample 

of 79.9 fb"^ recorded by the BABAR detector. We use the Cabibbo-favored D+ K'K'^tv'^ 
branching fraction as normalization in the measurements to reduce systematic uncertainties. The 
CP asymmetries obtained are AcpiK~ K+n^) = (1.4± 1.0(stat.) ±0.8(syst.)) x 10"^ Acp(07r*) = 
(0.2±1.5(stat.)±0.6(syst.)) x 10'^, and Acp{K*°K^) = (0.9 ± 1.7(stat.) ± 0.7(syst.)) x 10"^ The 

relative branching fraction ^^+'2^^-^+^+) a-lso measured and is found to be (10.7 ± 0.1(stat.) ± 

0.2(syst.)) X 10"^. 



PACS numbers: 11.30.Er, 13. 25. Ft, 14.40.Lb 



INTRODUCTION 



particle identification (PID) and tracking: 



Singly Cabibbo-suppressed (SCS) £>-nieson decays are 
predicted in the standard model (SM) to exhibit CP- 
violating charge asymmetries of the order of 10~^ 0. 
Direct CP violation in SCS decays could arise from the 
interference between tree-level (Fig. ^) and penguin 
(Fig. decay processes. Doubly Cabibbo-suppressed 
and Cabibbo-favored (CF) decays are expected to be CP 
invariant in the SM because they are dominated by a sin- 
gle weak amplitude. Measurements of CP asymmetries in 
SCS processes greater than 0(10"^) would be evidence 
of physics beyond the standard model 0. 

We define the CP asymmetry by 



Acp 



\A\' 



A 



(1) 



where A is the total decay amplitude for decays and 
A is the amplitude for the charge- conjugate decays. Acp 
is non-zero only if there are at least two different decay 
amplitudes with a CP-violating relative weak phase and 
a CP-conserving relative strong phase due to final-state 
interactions. Eq. |^ can be expressed as an asymmetry 
of branching fractions. We assume that the total decay 
rates for I?+ and D~ are equal (CPT invariance). As- 
suming further that CF decays are invariant under CP, 
we use branching fractions for CF decays as normaliza- 
tion factors to reduce experimental systematics due to 




*K+K--iT+) 



B(D- 



A, 



CP 



b{dT- 



*K+K-Tr-) 



B{D+^K + K--!r+) 
B{D+^K+K-7T+) 



B(D- 



*K+K-TV-) 



(2) 



(Throughout this paper we assume that the production 
of and _D+ mesons is charge symmetric.) 

We also measure the CP asymmetry in the resonant 



decays — > (pir^ and 



+ K*^K'^ , and determine 
the relative branching fraction ^y{d+^^^ -t+ ^+ ) ■ 



DETECTOR AND DATA SAMPLE 

This analysis is performed with a data sample recorded 
on and below the T{4S) resonance with the BABAR detec- 
tor at the PEP-II asymmetric-energy e+e^ storage rings 
at the Stanford Linear Accelerator Center. 

The BABAR detector is described in detail in Ref. 01. 
The silicon vertex tracker (SVT) and the 40-layer cylin- 
drical drift chamber (DCH) embedded in a 1.5-T solenoid 
measure the momenta and energy loss {dE/dx) of 
charged particles. A ring-imaging Cherenkov detector 
(DIRC) is used for charged-particle identification. Pho- 
tons are detected and electrons identified with a CsI(Tl) 
electromagnetic calorimeter (EMC). 

We split the 89.7 fb~^ data sample into a randomly 
selected subsample of 9.8 fb~^ to optimize the selection 
criteria and the remainder (a 79.9 fb"^ sample) for the 
final analysis. This procedure eliminates selection bias. 
We apply the same selection criteria to the CF and SCS 
modes whenever possible. We determine selection effi- 
ciencies from a sample (145 fb~^ equivalent) of Monte 
Carlo (MC) jj] generated e+e~ cc events. 



(a) 



(b) 



DATA ANALYSIS 



FIG. 1: Parton-level diagrams for D"*" — > K K'^n'* 
(a) a tree diagram, and (b) a penguin process. 



decays: -y^g reconstruct £)+ and Df Q decays by selecting 

events containing at least three charged tracks. Tracks 
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are required to have at least 12 measured DCH coordi- 
nates, a minimum transverse momentum of 0.1 GeV/c, 
and to originate within 1.5 cm in xy (transverse to the 
beam) and ±10 cm along the z-axis (along the e~ beam) 
of the nominal interaction point. Kaons are identified 
by a selection on the ratio of likelihood functions de- 
rived from dE/dx in the SVT and DCH, and from the 
Cherenkov angle and number of photons in the DIRC. 
Pions are identified as tracks that fail a loose kaon iden- 
tification criterion. The three charged tracks are further 
constrained to originate at a common vertex, the fit for 
which is accepted if the satisfies P(x^) > 1%. We 
reject Z?"*" and mesons from B decays, and thereby 
reduce backgrounds, by requiring that their momenta in 
the center-of-mass (CM) frame be above 2.4 GeV/c. 

In order to reduce the remaining combinatorial back- 
ground we consider likelihood ratios formed from the 
probability density functions (PDFs) of the following dis- 
criminating variables for the and D+ decays: CM 
momentum (pcm), vertex- fit probability with a beam- 
spot constraint (Pbs(x^))) ^ind the distance in the xy- 
plane from the interaction point to the D"*" or D'^ vertex 
[dxy]- The quantity i-Bs(x^) is the probability that the 
decay tracks form a vertex within the beam spot region. 
Most of the _D+ mesons decay outside this region, thus 
the probability P-esix^) is small for the Z?+ signal and is 
large for combinatorial background. Background distri- 
butions are taken from sidebands in the K~ K'^tt^ mass, 
while signal distributions are obtained from the signal 
regions, with the normalized sideband distributions sub- 
tracted. 

For decays, the rriKK-n signal band is defined as 
[1.840,1.896] GeV/c^ and the sideband mass regions as 
[1.805, 1.833] GeV/c2 and [1.903, 1.931] GeV/c^ [seeFig.El 
(a)]. Product likelihoods are constructed for the sig- 
nal, £sig = and the background, £bkg = 
rij '^bkg(^Oi where i runs over two or more of the vari- 
ables described. 

About 16% of the events have more than one _D+ me- 
son candidate. For such events the candidate with the 
highest likelihood ratio is selected. 

The sensitivity S//S.S, where S and AS" refer to the 
signal yield and its uncertainty, is optimized as a func- 
tion of the product likelihood ratio r = /Isig/'Cbkg formed 
using pcM and -Pbs(x^)! the optimal selection is found to 
be r > 4.3. This criterion is applied to both CF and SCS 
decays. When we use the analogous ratio ri obtained by 
including the PDF for d^y in £sig and £bkg, the sensitiv- 
ity is nearly as good. The results we find using ri provide 
a measure of systematic uncertainty. 

The subsamples for the decays Z?+ (pn^ and 
Z?"*" — > K*'^K'^ are selected by requiring that the invari- 
ant mass of the resonant decays be within 10 MeV/c^ and 
50 McV/c^ of the nominal (f> and K*^ masses, respec- 
tively |6|. In addition, the resonant signal samples are 
enhanced by a selection on the cosine of the helicity an- 



gle (cosOh)- In the (/>7r+ decay mode, the helicity 

angle is defined as the angle between the K~ and the tt"*" 
in the (f> rest frame. In the K*^K^ decay mode, 

the helicity angle is defined as the angle between the K~ 
and the K'^ in the K rest frame. Maximum sensitivity 
is obtained when |cos0_f/| > 0.2 and [cosS/fl > 0.3 for 
_D+ (Jyn^ and for D+ K*°K^, respectively. 

The CF — > K^K^n^ decays are selected by a pro- 
cedure identical to that for the SCS D+ K-K+it+ 
decays. We choose the signal rriKK-K region to be 
[1.944, 1.992] GeV/c^, while the sidebands are chosen 
to be [1.914,1.938] and [1.998,2.022] GeV/c^, respec- 
tively [see Fig. |5](b)]. In addition, contamination from 
K~Tr~^7r~^ decays is removed as follows: for all 
KKtt candidates, the kaon with the same charge as the 
pion is treated as a pion and then the Kttt: invariant 
mass is calculated. We observe a D'^ peak, indicating 
that part of the signal is composed of misidentified 
candidates. Events in the region 1.855 < mx-K-K l£ 
1.883 GeV/c^ are removed from the sample. Contam- 
ination from D*+ D°{-^ K-tt+ .K' K+)tt+ decays 
is removed by eliminating events for which m^-^+ > 
1.84 GeV/c^. Candidates for D+ K^tt+'k+ are elimi- 
nated if both Kit combinations satisfy the requirement. 
Partially reconstructed D*+ D°{^ K-ti+tt°)ti+ de- 
cays can also be misidentified as K~K^tt^ candidates if 
the 7r° is missed and the charged pion is misidentified as 
a kaon. Most of these decays are eliminated by assigning 
a pion mass to kaon tracks and removing candidates for 
which the mass difference {rux- tt+tt+ ~ 'T^_K'-Tr+) lies in 
the range [0.139,0.150] GeV/c^. 

Figure |21 shows the mass distributions obtained after 
all selection criteria are applied. The yields, listed in 
Tabled are computed by subtracting from the number of 
events in the signal region a scaled background estimate, 
obtained from the sideband mass region. 

TABLE I: Yields of background subtracted events, separately 
for each charge. 



Parent Charge + — 

K'K+n^ 21632 ± 228 20940 ± 226 
^ 07r± 5452 ± 87 5327 ± 86 

D^^K'^^K^ 5247 ±96 5113 ± 96 
Df K~K+Tv^ 23066 ± 217 22928 ± 214 



The efficiencies needed for the Acp calculation are ob- 
tained from a sample of MC generated cc events to which 
the same selection criteria are applied. The efficiencies 
for each decay mode are shown in Table Hll 

We obtain Aqp using Eq. (01 and replacing branch- 
ing fractions with efficiency-corrected yields. The results 
are shown in Table UlTl We also studied the CP asym- 
metry in 16 bins of the K~K~^7r~^ Dalitz plot 
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and found that the asymmetry is consistent with being 
constant (with a probabihty of 51%) and zero. 

We use the CF sample of £)+ K^tt^tt^ de- 
cays, obtained using selection criteria identical to the 
SCS case, to determine the relative branching fraction 

^^^)+^^-^4-^+^ as follows. The CF and SCS Dalitz plots 
are first divided into equally populated bins (16 bins for 
the SCS mode, 64 for the CF mode). Next, the signal and 
normalization yields and efhciencies are calculated bin- 
by-bin. The efficiency-corrected yields are then summed 
and divided to obtain the ratio. Figure O shows the mass 
distribution in the CF — > K^tt^tt^ mode, for which 
the average efficiency is 10.03 ± 0.01(stat.)%. We obtain 
a relative branching fraction of (10.7±0.1(stat.)) x 10"^. 

TABLE II: Efficiencies for positively (e"*") and negatively (e~) 
charged D and Ds meson decays. Efficiencies are in percent. 
The stated uncertainties are due to MC statistics only. 



Decay e 
R-R+n^ 8.20±0.04 8.26±0.04 
^ <?!)7r± 7.67±0.07 7.63±0.07 

D± ^ K*'^K^ 5.88±0.07 5.90±0.07 
Df K~K+-!v^ 3.77±0.02 3.79±0.02 



asymmetries in £)+ and Df decays. We evaluate an un- 
certainty for the background subtraction by changing the 
widths of the sideband mass regions. The uncertainty is 
taken to be the difference in the central values of Acp- 
The uncertainties in the likelihood-ratio technique are es- 
timated with two variants: (i) tightening the likelihood 
ratio to produce a 10% change in the yields, and (ii) us- 
ing the likelihood ratio ri (described above) in place of r. 
The systematic uncertainty is chosen to be the larger of 
the two changes. Table HVl summarizes these systematic 
uncertainties for the observed CP asymmetries. 



TABLE IV: Systematic uncertainties for the CP asymme- 
tries. 



Source 




(^7r± 






Acp [10-2] 


Acp [10-2] 


Acp [10-2] 


MC simulation 


0.06 


0.06 


0.06 


Background estimate 


0.63 


0.32 


0.49 


Selection criteria 


0.46 


0.54 


0.54 


Total 


0.78 


0.63 


0.73 



TABLE III: Results of the CP-asymmetry measurements, 
Acp- Also listed are the values for A^^p, the asymmetry 
computed without the normalization mode. 



Decay Acp [10~2] A'-^^ 

K-R+Tv^ +1.36 ±1.01 +2.07 ±0.84 

^± +0.24 ± 1.45 +0.94 ± 1.33 

K*°K^ +0.88 ± 1.67 +1.58 ± 1.57 



SYSTEMATIC UNCERTAINTIES AND 
CROSS-CHECKS 

The only difference between the final states from Df 
and D^ decays considered here is a slightly harder mo- 
mentum spectrum for the Df decay products. In turn, 
these small differences are corrected for by the efficien- 
cies which come from MC. Any charge asymmetry in the 
detection of pions thus cancels when Df K~K~^tt^ 
decays are used as normalization, as in Eq. 0. We esti- 
mate the systematic uncertainty on the CP asymmetries 
by combining estimates of the contributions from various 
identified sources listed in Table HVl 

The uncertainty due to small differences in momentum 
spectra of tt, K from £>+ and Df decays, 0.06%, is es- 
timated as three times the maximum difference in tt, 



We performed two cross-checks on our measurement 
of Acp- First, we calculated an alternative measure of 
CP asymmetry without using Df K~K~^Tr~^ decays 
as normalization, which we labeled ^[?], in Table Hn] 
We find its values to be consistent with our measure- 
ments of Ac p ■ Second, we measured the CP asymmetry 
for a control sample: the CF decays Df K^K^tt'^ 
(non-resonant as well as resonant). This asymmetry is 
expected to be zero within the SM. In Df K~K^tt^ 
decays, both the Df and the Dj decay to two oppositely 
charged kaons and only the pion charge differs in parti- 
cle and antiparticle decays. Thus, any detector-induced 
asymmetry would arise only from a charge asymmetry in 
pion tracking and is expected to be very small. Indeed 
the measured value is (+0.6 ± 0.8) x 10~^. 

As a final cross-check, the CP asymmetry has also been 
studied as a function of the D+ laboratory momentum, as 
well as by the run period. No significant dependence on 
momentum or detector operation conditions is observed. 

A summary of the systematic uncertainties for the rel- 
ative branching fraction -"^+ ^+ ^ is given in Ta- 
ble The fractional uncertainty due to PID and track- 
ing has been estimated as 2.1%, computed as the sum 
in quadrature of 1.1% for PID and 1.8% for tracking 0|- 
The PID uncertainty is estimated from a comparison of 
PID efficiencies in data and MC. The tracking uncer- 
tainty, which is the uncertainty on the K/ir efficiency 
ratio, is conservatively estimated as three times its value 
obtained using MC. 
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FIG. 2: KKn mass distributions for positively charged (left) 
and negatively charged (right) D and Ds candidates for events 
satisfying the requirement r > 4.3. Figures (a) and (b) are 
for all KKn candidates, while (c) and (d) are for <j!>7r candi- 
dates, and (e) and (f) for K K candidates. Signal (yellow 
or light shaded) and sidebands (red or darker shaded) regions 
are shown for and Df decays in (a) and (b), respectively. 
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TABLE V; Systematic uncertainties for the relative branch- 
ing fraction. 



Source 


Uncertainty [10 


PID -f- tracking 


0.22 


Background estimate 


0.05 


Selection criteria 


0.02 


Total 


0.23 



TABLE VL Results of the CP asymmetry (Acp) measure- 
ments for decays. 



Decay 


Acp [10-"] 


K-K+n^ 


-f 1.4 ± 1.0(stat.) ± 0.8(syst.) 
+0.2 ± 1.5(stat.) ± 0.6(syst.) 
+0.9 + 1.7(stat.) + 0.7(syst.) 



onant modes having significantly smaller uncertainties. 

Further, we obtain a branching fraction for 
K'K^TT'^ decays relative to that for _D+ K^tt^tt^ 
decays of (10.7 ± 0.1(stat.) ± 0.2(syst.)) x 10"^^ This 
result is a significant improvement over previous mea- 
surements 113]. 
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INFN (Italy), EDM (The Netherlands), NER (Norway), 
MIST (Russia), and PPARC (United Kingdom). Indi- 
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